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ABSTRACT 

We report on the highest to date signal-to-noise-ratio X-ray spectrum of the lu- 
minous quasar PDS 456, as obtained during two XMM-Newton orbits in September 
2007. The present spectrum is considerably different from several previous X-ray 
spectra recorded for PDS 456 since 1998. The ultra-high-velocity outflow seen as 
recently as February 2007 is not detected in absorption. Conversely, a significant re- 
flection component is detected (Ax^ = 313 compared to a simple absorbed power law). 
The reflection model suggests the reflecting medium may be outflowing at a velocity 
v/c = -0.06 ± 0.02 (Ax^ = 28 compared to v/c - 0). The present spectrum is ana- 
lyzed in the context of the previous ones in an attempt to understand all spectra within 
the framework of a single model. We examine whether an outflow with variable par- 
tial covering of the X-ray source along the line of sight that also reflects the source 
from other lines of sight can explain the dramatic variations in the broad-band spectral 
curvature of PDS 456. It is established that absorption plays a major role in shaping 
the spectrum of other epochs, while the 2007 XMM-Newton spectrum is dominated by 
reflection, and the coverage of the source by the putative outflow is small (< 20%). 

Subject headings: quasars: individual: PDS 456 - X-rays: galaxies - techniques: 
spectroscopic 
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1. Introduction 



PDS 456 was discovered as a radi o-quiet luminous c 



survey (PDS) of bright IRAS sources ([Torres et al. 



uasar in the course of the Pico dos Dias 



19971). The optical spectrum of PDS 456 



features broad emission line s and a redenning-corrected luminosity 30% higher than that of 



3C 273 ([Simpson et al 



19991) ■ The high luminosity of PDS 456 provides a low limit of 10'^ Mc 



on the black hole mass, just from the Eddingt qn limit. Rapid X-ra y variability of PDS 456 was 



first noticed with the RXTE satellite in 1998 by 



Reeves et al.. (i2000t) . who observed the 2-10 keV 



flux double in 17 ks, which is rather uncommon for such a luminous quasar and indicative of an 
unusually compact source of a few gravitational radii. 



A few years later, 



Reeves et al. 



(120031) reported an ultra high velocity outflow from PDS 456 



based on blue-shifted Fe K-shell absorption features detected in a 40 ks XMM-Newton observation 



from 2001. Extrernely b 



(e.g., 



Chartas et al 



2002, 



ue-shi f ted Fe-K absorption has been reported for other quasars as well 



2003 



Pounds et al. 



20031), Obviously, for a confirmed detection of 



an outflow and for conclusive velocity measurements, it is most desirable to observe as many 
resolved lines as possible, all shifted by the same velocity. This is most difficult, since in 
many objects the absorption is detected primarily from the most highly ionized species of Fe, 
whose lines occur at energies where the sensitivity of current detectors drops sharply, and in 
ionization regimes that preclude absorption by all other species that are too (fully) ionized to 
produce lines. The difficulty is then augmented by the rapid sp ectral variability of these outflows 



Kaspi & Behar 



20051: 



Chartas et al. 



20071 : 



Reeves et al. 



20081) . which imposes fundamental limits 



on the signal-to-noise-ratio (S/N) of spectra and deems data stacking from separate observations 
almost useless. All these have raised questions regarding t he statistical robustness o f some, but 



probably not all, of the blue- shifted absorption detections (|Vaughan & Uttley 



2008b . The only 



way to substantiate and further study these high velocity outflows is with long exposures and high 
S/N spectra, such as the present one. 
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PDS 456 stands out as one of the least controversial ultra h igh velocity outflows. Not only 



2000 



2003 



2009h . 



does it have multiple Fe-K absorption detections over the years (iReeves et alj 
but a 2001 XMM-Newton observation revealed an absorption feature in the soft X-ray band that was 
measured at high-resolution with the Reflection Grating Spectrometers (RGS). Indeed, the RGS 
spectrum sho ws a deep and broa d absorption trough between 13-14 A (in the observed frame), 
interpreted by .Reeves et al..(i2003.) as due to a ~50,000 km s"^ outflow that also produces the Fe-K 
absorption edge. Blue-s hifted Ly a absorptio n indicates an appreciably slower outflow moving at 



14,000 - 24,000 km s'^ dO'Brien et al 



120051) . Despite its moderate redshift of z = 0.184, PDS 456 
is also the brightest among the targets suspected of harboring an ultra high velocity outflow. With 
all this in mind, we obtained a deep XMM-Newton exposure of PDS 456, which is reported in this 
paper. XMM-Newton has the largest efi'ective area of all current X-ray telescopes and thus provides 
the best chance to observe several absorption lines from the same outflow component with the 
gratings. 

In the meantime, a Suzaku observation of PDS 456 provided another high-quality (though 
low-resolution) spectrum of PDS 456 in which Fe-K line shi fts of 100,000 km s ^ were detected 



at high S/N, but with no obvious soft X-ray line absorption (IReeves et al. 



2009). This adds to 



the mounting evidence of blue-shifted absorption detected in PDS 456 above 8 keV. An exciting 
discovery in the Suzaku observation is the significant emission at high energies above 15 keV, 
much in excess of the extrapolated soft X-ray spectrum. This is yet another manifestation of the 
unique spectral behavior exhibited by this intriguing quasar. In the present paper, after analyzing 
the 2007 observations and studying the short term X-ray variability within those two XMM-Newton 
orbits in Section[3l in Section|4]we also study the spectral variability of PDS 456 between different 
epochs, and on time scales of months and years, over the past decade. 
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2. Observations 

2.1. Data Reduction and Light Curves 

In Table [H we list the log of the 2007 XMM-Newton observations of PDS 456 as well as the 
earlier XMM-Newton observation from 2001, and then the other observations of PDS 456 that are 
used later in this paper. XMM-Newton data were reduced using the Science Analysis Software (SAS 
- ver. 7.1.0) in the standard processing chains as described in the ABC Guide to XMM-Newton Data 
Analysi^, and using the most updated calibration files for 2009, June 15. Both single and double 
events were retained and no pile up was detected. 1.5 ks were excluded from the orbit- 1421 data 
set due to high background. 

The light curve of the source during orbits 1421 and 1422 is shown in Fig.[TJ The source 



flux diminishes gradually during the second orbit to approximately half its initial flux 



Reeves et al. 



evel. The 



(120021) . and 



typical time for variability in Fig.[T]is ~30ks, similar to that reported by 
which corresponds to the light crossing time of only a few gravitational radii. It should be noted 
that despite the significant X-ray variability of PDS 456, it is rather steady in longer wavebands. 
The optical monitor (OM) on board XMM-Newton recorded a flux range of (2.4 - 2.5) xlO 
ergs"^ cm"2 (1.10 - 1.14 cts s'^ during orbit 1421) with the UVM2 filter around 2300 A and 
a range of only (5.75 - 5.85) xlO~'^ erg s~' cm"^ (12.1 - 12.3 cts s"' during orbit 1422) with 
the UVWl filter around 2900 A. These can be further compared with the archival HST/STIS 
measured flux of 3 xlO"'^ erg s"^ cm'^ at 2300 A and 9.5 xlO"^^ erg s*^ cm'^ at 2900 A on 2000, 
May 14. All fluxes are quoted before correcting for reddening. Optical monitoring in the V band 
carried out at the Wise Observatory every two weeks over half a year reveals an increase in flux 
from 2007 April to July of only 0.05 mag, and then a decrease in flux from 2007 July to October 
of 0.07 mag. Thus, while the X-ray flux varies by a few tens of percent over a period of a day. 



^See http://heasarc.gsfc.nasa.gov/docs/xmm/abcl 
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the UV flux is constant to within a few percent. At much longer time scales, perhaps even years 
but not well monitored, the UV flux can also vary by a factor of two. The optical flux varies by 
only a few percent over several months. This is not unusual, as quasar monitoring often shows 
an increase of variab ility with frequency from optical, through UV, and up to the X-ray band 



( Edelson et al 



1996 ). 



2.2. Background 

In order to obtain the X-ray spectra, we extracted the source photons from the EPIC PN 
and MOS cameras in a circular region of radius 40" and grouped the data to bins with minimum 
25 counts per bin. We use the standard SAS background files for each observation, where the 
background spectrum is extracted from an annulus region centered around the source of 1' to 
3', from which 1' circular regions around field point sources are excluded. The background 
for EPIC PN requires special caution as it can have strong instrumental, narrow emission lines 
at energies around 1 .5 and 8 keV, which could introduce by subtraction artificial absorption 
lines. The background during the 2007 observations was relatively high, but we verified that the 
background spectra we use have no such lines. 

We measured the background in various off'-source locations on the EPIC PN and MOS 
CCDs. A sample background spectrum is shown in Fig. [2l Indeed, instrumental emission lines 
appear on th e detector with varyin g strengths as has been documented in the XMM-Newton Users' 



Handbook 



Freyberg et al 



(|2004|) investigated the intrinsic instrument background in the EPIC 
detectors and found that it is composed of a continuum and X-ray fluorescence lines, most 
prominent are those at 1.5, 7.7, 8.0, 8.6, and 17.4 keV, which are due to Ka line emission of Al, 
Ni, Cu, Zn, and Mo, respectively. The lines show strong spatial inhomogeneities correlated with 



^http : //heas arc . gsf c .nasa.gov/docs/xmm/uhb/node3 5 .html 
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the structures of the electronics board mounted below the EPIC PN CCDs. iFreyberg et alj (|2004|) 
also found that the very central region of the camera where typically the primary source is located 
is free of this instrumental background. Therefore, our use of line free background spectra is 
justified. 

We also examined the background in the archival XMM-Newton 2001 observation. For this 
observation as well, we use the standard 1' to 3' annulus. The background in the 2001 observation 
is much lower and its fluorescence lines much less prominent. More precisely, at the edge of 
the CCD the fluor escence lines are stro ng, but they get much weaker close to the center of the 



CCD, as expected (IFrevberg et al 



20041). Thus, the spectral turnover detected around 8 keV in the 



2001 spectrum of PDS 456 (IReeves et al. 
background. 



2003|) is secure and not significantly influenced by the 



3. 2007 XMM-Newton Campaign 

The September 2007 XMM-Newton observations provide the highest S/N spectra of PDS 456 
to date. EPIC PN, for example, collected an estimated 437,600 source photons, compared to 
121,700 photons during the 2001 observation. The spectra of both EPIC cameras and the RGS are 
of sufliciently high statistical quality to be confronted with a variety of models. Below 0.4 keV, we 
find discrepancies of 20% between the EPIC PN and MOS data with respect to a simple power-law. 
Hence, the present modeling is carried out in the 0.4-12.0 keV range and highlights mostly the 
more sensitive EPIC PN. These spectra are presented separately for the XMM-Newton orbits 1421 
and 1422 in Fig. [31 along with fits to increasingly complex models. Although the separation of 
the data into the two XMM-Newton orbits is arbitrary, it provides an insight into the short-term 
variability and a reality check for persisting spectral features. It is noted that the EPIC PN-fitted 
models of each orbit also fit the MOS spectra very well. For presentation purposes alone, we have 
omitted the MOS data from the plots. 
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We first fit the spectra with a simple power law absorbed by the nom inal neutral H Gala ctic 



Wilms et al 



(l2000h 



column density of Nh = 2.0 x 10^^ cm"^, using the absorption model of 
We find that the photon slope Y slightly increases by 0.05 between the two orbits. This model 
alone is clearly inadequate, as seen by the large residuals plotted in the top panel of Fig. [31 
and since it yields ^^d-O-f = 1-9. The residuals below 1 keV suggest the neutral column is 
underestimated. Indeed, the best-fit neutral column density is Nh = (2.67 + 0.03) x 10^^ cm"^, 
which considerably improves the fit to;^Vd.o.f =1.19 (middle panel of Fig. [3]). Separate fits to 
the RGS and MOS spectra yield respective column densities of (2.64 + 0.10) x 10^' cm~^ and 
(2.65 + 0.05) x 10^' cm"^, all consistent and robustly above the Galactic column. The redshift 
of PDS 456 (z = 0.184) allows us to test whether the excess column occurs in the quasar, by 
invoking a redshifted absorption component in the model. However, the fit decisively ascribes the 
entire neutral column to the z = component and no column to the redshifted one. Moreover, if 
the redshift of the additional component is left to float, the fit requires that it converges to z = 0. 
The agreement between EPIC PN, MOS, and RGS imply that dete ctor contamination can not be 



the reason for the excess absorption (see also 



Kirsch et al. 



20051). We tested whether perhaps 



a partially covering absorber can explain the excess low-energy absorption. However, the best 
fit model of a partial covering neutral absorber yields a covering fraction of 100% (and with no 
improvement to;(fVd.o.f = 1.19). There is also no indication of an ionized absorber that could be 
responsible for the excess absorption, as no absorption lines are detected anywhere in the spectra. 
Indeed, the best fit for a partial covering ionized absorber yields a column density that tends to 
zero, and again no improvement to;(fVd.o.f = 1.19. 

To make a long story short, the data clearly indicate excess neutral absorption that is 
consistent with z = 0. Although in principle it is possible that by coincidence this absorption is 
due to gas in PDS 456 that happens to be outflowing by exactly the cosmological recession velocity 
of the quasar, we reject this notion as highly unlikely. Moreover, since all indications point to 
this absorber being neutral, we do not associate it with the highly ionized outflows observed in 
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PDS 456. Similar variable Galactic absorption has been detected towar ds NRAO 140, wh ere it 



19951) . Does 



could be explained by motion in the Perseus molecular cloud complex (|Tumer et al. 
the moderately low Galactic latitude of PDS 456 of +11° suggest a similar explanation in the 
present case? We searched for an y molecular cloud s towards PDS 456, both in the IRAS archive. 



and in the Columbia CO survey (IDame et al 



1987h . but found none. The origin, therefore, of 



the excess neutral absorption and its discrepancy with previous measurements remains somewhat 
unclear. 

In the middle panel of Fig. [3l we show the spectral residuals that remain after incorporating 
the higher column density in the model. It can be seen that these residuals are fairly similar in both 
orbits (but keep in mind the slightly different power law slopes). The residuals below 1 keV could 
be due to the incomplete absorption model around the 0-K (~ 0.53 keV) and Fe-L (~ 0.71 keV) 
photoelectric edges. At high energies, there is excess emission around the Fe-K region between 5 
and 7 keV in the observed frame. Note that much of this emission is blueshifted with respect to 
the quasar frame, given that the 6.4 - 7.0 keV range of Fe-K lines corresponds in the observed 
frame (Fig. [3]) to the 5.4-5.9 keV energy range. The Fe-K excess emission and its high statistical 
significance is further demonstrated in the;^'^ residuals plotted in Fig. HI The second data set 
(orbit 1422) has a narrow component, but both data sets feature large broad residuals resembling 
the electron-scattered blue wings typical of ionized (Compton) reflection. Also important is the 
over all curvature of the model between 1 and 5 keV, where no narrow features occur, that is not 
properly reproduced by an absorbed power law. Together with the Fe-K residuals this strongly 
suggests a reflection spectral com ponent. Indeed, adding a constant density reflection component 



(reflion, by 



Ross & Fabian 



120051) notably and significantly improves the agreement of the model 
with the data (Ax^ = 313, and;i^^/d.o.f = 1.05). The broad-band spectral curvature between 1 
and 5 keV is also much improved when reflection is included in the model, as can be seen by 
comparing the middle and bottom panels of Fig. [31 The small remaining residuals are of order of 
10 - 15 %. The best-fit parameter values for the two combined orbits are given in Table [2l When 
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fitting the two data sets separately, we find that the reflection component fits the second data set 
(orbit 1422) particularly well, yielding ;^Vd-o.f = 0.99. The fit to the first spectrum (orbit 1421) is 
somewhat poorer with;^Vd.o.f = 1.10. Models fitted separately to each spectrum yield parameters 
that are similar to those of Table [2] (log ^ « 4,AFe « 3,z ~ 0.1), but are not as tightly constrained. 

We tested alternative models to reflection, including a plain Gaussian emission line and a 
gravitationally distorted emission line. Both give unphysical results as could have been expected 
from the spectral shape of the residuals and the combined narrow and broad nature of the residuals 
in the middle panel of Fig. [3l A Gaussian fit forces an enormously broad feature with cr = 
3 keV, and centered at an unreasonably low energy of « 4 keV in the quasar frame, essentially 
mimicking the underlying broad reflection continuum. A relativistically distorted line profile 



(Laoi 



19911) also yields unreasonable physical parameters, such as a disk inclination of almost 
90° and an emissivity that increases with radius. Moreover, partial covering neutral, and partial 
covering ionized absorption models were also tried in place of reflection, but could not eliminate 
the residuals observed in the middle panel of Fig. [31 In summary, we prefer the reflection model, 
which is most appropriate for describing both the broad-band and the Fe-K spectrum of PDS 456 
during the 2007 XMM-Newton observations. 

The fit of the reflection model to the combined 2007 spectrum constrains the frame of 
reference of the reflection component rather tightly. This comes from the position of the reflected 
Fe line (for log^ ~ 4 primarily H-like at 6.97 keV) and even more from the electron-scattered 
blue wings of the lin e that give it a slightly broadened and skewed profile (see, e.g.. Fig. 3 in 



Ross & Fabian 



2005 



and also Fig.|9]below). Interestingly, the best-fit redshift for the reflection 
model is indicative of a frame of reference that is distinct from both the quasar and the observer, 
namely z = 0.12 + 0.02 to be compared with z = 0.184 for PDS 456. Fixing the reflection redshift 
to that of the quasar and refitting results in an inferior fit (A;^f^ = 28 resulting in;^f^/d.o.f. = 
1.07). The intermediate reference frame can be understood as reflection off of fast moving gas. 
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perhaps the outflows detected in previous observations? Since both the frame of reference and 
the ionization state of the reflecting medium bear on the position of its Fe-K line, it is important 
to investigate the connection between the two in the fit. In Fig. [51 we plot the combined statistical 
confidence contours for the redshift and the ionization parameter. The plot demonstrates that the 
kinematical shift (i.e., outflow) of the reflector is statistically robust (> 3cr) against the ionization 
uncertainties within the framework of the present model. However, this result is model dependent 
and a realistic outflow is surely more complicated than the slab reflection used here. The regime 
of ionization parameter beyond ^ = 10,000 (cgs units) is inaccessible with the current reflection 
models, and therefore is not plotted in Fig. [51 Moreover, the constraints on each XMM-Newton data 
set when fitted separately are insufficient to produce similar contour plots, only the combined fit is 
well constrained. On the other hand, much higher ionization > 10,000 in cgs units) can not shift 
the line position to beyond its current H-like position, which suggests that even higher ionization 
would perhaps not invalidate the outflowing nature of the reflector. There is also no obvious 
indication by extrapolation of Fig. [51 that higher ionization compromises the redshift. The putative 
radial velocity of the reflector, which would be moving more or less parallel to but outside the line 
of sight is Vout/c ^ (0.12 + 0.02) -0.184 = -0.06 + 0.02, or Vout = -18,000 + 6, 000 km s^^ In 
the next section, we show that the moving reflection component may be a persistent feature of the 
otherwise variable X-ray source in PDS 456. 

Furthermore, if the observed reflection component is due to reflection of the observed 
primary power law, it is possible to estimate the solid angle occupied by the reflector as viewed 
from the primary source. In a simplified geometry with perfect reflection (i.e., no flux losses on 
the reflector), and neglecting time lag, the flux ratio of the reflection component to the primary 
power law gives a rough idea of this fractional sky coverage Qj-^-fi/ln. The n o-losses approx i matio n 



Ross & Fabian 



mm 



is reasonable for electron reflection in a nearly fully ionized reflector. The 
model assumes uniform illumination of the reflector and a plane parallel geometry. It computes 
the reflection integrated over all angles from 1 eV to 1 MeV by a power law source ranging from 
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1 eV to 0.3 MeV. For the present models, the flux ratio in the observed region from 0.4 - 12.0 keV 
gives essentially the same result as that in the full range. For the best fit model, this flux ratio 
implies Qrefl/27r « 0.2. The best-fit reflected flux in the present models happens to be similar in 
both orbits, despite the decrease by approximately 30% of the mean primary power law flux. This 
indicates the reflector is likely extended from the source. 

Given the multi-component nature of the PDS 456 X-ray spectrum, and given the variability 
of shifted emission and absorption on short and long time scales, as well as the low Galactic 
latitude (-1-11°) of PDS 456, it is worthwhile to examine the viability of a previously unknown 
Galactic source along the line of sight. An X-ray binary (XRB) could be the primary candidate for 
such a source. In order to test such a conjecture, we inspected the X-ray Chandra image, the UV 
XMM-Newton optical monitor (OM) image, and the 2MASS K-band image of PDS 456. They are 
all perfectly consistent with point sources and are co-aligned to within a fraction of an arcsecond. 



Grimm et al 



(l2002h for galactic XRBs down to 



By extrapolating the log N - log S distribution of |j 
a flux, say, of / = 10"'^ erg s~' cm"^ (about 10% of the current flux of PDS 456, but two orders of 



Grimm et al. 



20021), one can use the expression 



magnitude fainter than the completeness limit of ll 
A^(> S) = 72(//3.2 X 10"^° erg s"' cm"2)-o.4i estimate the number of XRBs in the sky that 
are brighter than 10"'^ erg s ' cm"^ to be not more than 800. Consequently, the probability for a 
random XRB within the solid angle of 0.5"^ is < 4 x 10^'°; negligibly small. In other words, it 
seems unlikely that a foreground XRB is present along the line of sight to PDS 456. 

In Section |4l we explore partial covering models by outflows as a way to explain on an 
equal footing the entire slew of X-ray spectra obtained for PDS 456 over the past decade. In the 
2007 EPIC PN spectra shown in Fig. [3] there is no obvious absorption by an outflow, and thus the 
absorption models actually do not improve the fit. The reflection component, on the other hand, 
seems to persist in all epochs. 
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3.1. Outflow? 



We searched the present 2007 data set for the fast outflows detected by 



Reeves et al. 



(12003 . 



20091). The high-energy region of the EPIC PN spectrum during orbit 1421 is shown in Fig. |4l 
No sign of absorption can be detect ed whatsoever; Nei ther the lines that appear in the Suzaku 



spectrum at 7.67 keV and 8.13 keV (IReeves et al. 



XMM-Newton observation (Reeves et al 



20091), nor the edge that appeared in the 2001 



20031) are seen in this data set. We further inspected the 



RGS spectra from both XMM-Newton orbits. Fig. [6] presents the ratio of the RGS spectra to the 
model used above to fit the EPIC PN spectra, with no adjustments whatsoever. Since this model 
includes no discrete spectroscopic features (except perhaps the neutral O absorption at « 23.5 A), 
we would expect absorption (or emission) lines or edges due to the outflow to show up in the RGS 
data. We show in Fig. [6] each RGS (1 and 2) and the data of each orbit separately in an attempt to 
identify outstanding features. The statistics of the data presented is rather high, ranging from over 
5,000 to over 10,000 source photons in each panel, which are distributed over « 170 bins. In fact, 
this is the best high-resolution X-ray spectrum of PDS 456 to date. A closer look into the 13-20 A 
spectral region, where the RGS is most sensitive, is presented in Fig.|71 We are unable to identify 
any recurring lines that appear in all of the panels of Figs.[6]or|71 At the very least, we would 
expect genuine absorption lines to be present in both RGSs during a single observation (panels 1 
and 2, or 3 and 4), but we find none. We conclude that there is no unambiguous sign in September 
2007 of the outflows detected in the two previous observations. The intermittent nature of the 
outflow has immediate implications on the overall AGN energy budget as affected by such flows. 



4. Spectral Variability 

The appearance and disappearance of the outflow in PDS 456 within several months is 
only one aspect of the full range of changes taking place in this intriguing quasar. Since it was 
discovered, PDS 456 was observed by almost all available X-ray telescopes. Therefore, one can 
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expect to obtain a more comprehensive picture of this source by studying its X-ray spectral history. 
In addition to the current data sets from the two XMM-Newton orbits in 2007, there are five archival 
X-ray observations of PDS 456. These are available from XMM-Newton (February 2001), Suzaku 
(February 2007), Chandra (May 2003), ASCA (March 1998), and RXTE (Febru ary-March 2001). 



With the ex ception of Chandra, all observations were reported in the literature (|Reeves et al.ll2000 



2003L 



2009|) . For the details about these data and how they were reduced we refer the readers to 
those references. The Chandra HETGS observation was not published because it caught PDS 456 
in a very low state. The data are reduced here using CI AO version 3.2.1 and CALDB 3.0.1 and the 
standard CIAO threads. Spectral fits were carried out simultaneously for both HETGS instruments 
(MEG and HEG), although we show only the MEG spectra. In Fig. [H we plot the ratio of the 
different PDS 456 spectra to a simple Galactically absorbed {Nh = 2 x 10^^ cm"^) power-law with 
photon index F = 2. From the 2007 XMM-Newton campaign, only the spectrum of orbit 1421 is 
presented, to minimize crowding of the plot. The total flux variability is of the order of a factor of 
2-3. This may not be expected from a quasar as bright as PDS 456 {Lx ~ 10"^' erg However, 



PDS 456 has already been observed varying by a factor of a few over less than a day (|Reeves et al. 



200Qh 



From Fig. [81 it can be seen that spectral structure below 1 keV varies both overall and in its 
narrow features. The Fe-K region beyond 5 keV also changes drastically. What appears as an 
absorption edge, and can be interpreted as an outflow, varies both in position and in depth. The 
significance of emission features, possibly identified as due to reflection, changes as well. Between 
1 keV and 5 keV, the spectral curvature changes in what seems to be an intricate combination 
of emission and absorption variability. This observed complexity can not be characterized by a 
single spectral component or one obvious spectral parameter, such as power-law slope, Fe-K line 
flux, or absorbing column density. What on the face of it seem to be common features in several 
spectra, such as an Fe-K absorption edge (or emission line?) beyond 6 keV, or an Fe-L emission 
hump below 1 keV, upon closer inspection in fact appear to occur at slightly diff'erent energies in 
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each spectrum. 



In the next section, we attempt to use one possible type of models to fit all of the spectra 
of PDS 456. For this exercise we fit separately six epochs with high quality CCD data at an 
energy resolution of x 100 eV, or better, namely those of EPIC PN (two XMM-Newton orbits in 
September 2007), EPIC PN (2001), Suzaku XIS & 3 (February 2007), Chandra HETGS (2003, 
MEG & HEG) rebinned to CCD resolution due to low S/N, and A5CA SIS (1998). We leave out 
the RXTE spectrum due to its lower resolution. Data from MOS on XMM-Newton and XIS 1 on 



Suzaku are less sensitive at hi gh energies, 
XIS & 3, respectively (see iReeves et al. 



DUt their data are consistent with those of EPIC PN and 



2003 



2009). The six data sets used here are listed in 



Table [TJ The modeling strategy is described in the next section. 



4.1. Partial Covering Approach 



For the model, we invoke a simple continuum power law source, partially covered by a 
(possibly ionized) absorber and absorbed locally in our Galaxy. For the absorb ing outflow, the 



200 1|) available in Xspec 



analytical warmabs models version 2.11n8 of XSTAR ( Kallman & Bautista 

w 

used A reflection component is also included as in Section [31 to account simultaneously 



are 



for Fe-K lines and for the soft excess that is observed in some data sets and apparently varies in 
its discrete features (Fig. [H). In Xspec code, the model used is TBabs {zpowerlw + zpowerlw * 
warmabs * warmabs + reflion). That partial cov ering of the central so urce can provide useful 



insights into the AGN geometry was realized by 



Reichert et al 



(|I986h . These types of models 



provide a wide range of broad band spectral curvature that may (or may not) be able to produce 
the vastly different spectra seen in Fig.[8l The model discussed above for the 2007 XMM-Newton 
spectra (Section [3]) where no absorption is detected qualifies as the degenerate extreme case of 



^http://heasarc. gsfc.nasa.gov/docs/xanadu/xspec/models/xstar.html 
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zero coverage of the source. Realistic partial covering models also include narrow absorption 
lines, which impose additional constraints complementary to the broad-band curvature of the 
spectrum. However, when;^f^ minimization techniques are used to select models, as in the present 
work, the continuum statistics overwhelms that of the discrete features. Thus, the narrow features 
are much harder to constrain with the current CCD spectra. S imilar models of pa rtial covering 
by an ionized outflow have been tried w ith relative success by .Reeves et al.l (l2003i) for the 2001 



observation of PDS 456, and recently by 



Miller et al. 



(l2008[) for MCG -6-30-15. 



The modest goal of this section is to test the conjecture that perhaps it is the varying fractional 
coverage of the source by the outflow that drives the spectral variability, and whether these models 
can provide an acceptable overall explanation to all of the six spectra of PDS 456. This type of 
partial coverage models is complex and invokes by far more parameters than can be constrained 
by the data. Consequently, it is beneficial to fix a-priori as many parameters as possible. The 
illuminating spectrum to be absorbed and reflected is the aforementioned power law with its 
slope fixed to F = 2.25, which is the best fit value for the high S/N and relatively featureless 
spectra obtained with EPIC PN (orbit 1422) and with Suzaku. Two ionization components are 
included for the absorber. These components represent high and low ionizat ion in an absorber 



that more realistic ally could have a continuous distribution of charge states (|Steenbrugge et al 



2003 



Behai 



20091) . and are therefore not well constrained by low-resolution data. The fractional 
coverage of the power law source by the outflow is a fitted parameter. The outflow velocity is 
fitted in each spectrum and is found to somewhat vary from one observation to the other. The 
velocity broadening is fixed at Vturb = 2,500 km s~^ which fits well the only grating measurement 



absorption lines in the Suzaku spectrum (IReeves et al. 



of absorption lines in PDS 456 (see Sec tion 14. 3lbelow) and is consistent with the width of the 



2009) 



Since the reflection component in the 2007 XMM-Newton spectra appeared to be moving 
(Section [3]) and inconsistent with the rest frame of PDS 456, it may be moving with the absorbing 
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medium. Consequently, we tie the outflow velocities of the reflection and absorption components, 
which would be a good approximation as long as the opening angle of the putative reflecting 
outflow is not too large. Abundances are set to their solar values, except for the Fe abundance of 
the reflector that is set at three times solar as required by the fits in Section [31 No independent 
Fe-K absorption or emission feature beyond the outflow and reflection features is added. Also 
for simplicity, no relativistic blurring is invoked, since it tends to make the models even more 
degenerate than they already are. Despite all of these impositions on the model, it is still fairly 
degenerate and the fit can not reliably constrain all of its free parameters. To that end, the Suzaku 
spectrum can be singled out with its detected absorption lines from tw o kinematic compon ents 



2009 ). 



and the ionization parameter that can be more accurately determined (|Reeves et al. 
Consequently, in our model for the Suzaku spectrum, the absorber (v/c ~ -0.3) and reflector 
(v/c ~ -0.12) are allowed to have different velocities. 



4.2. Fit Results 

Using the partial covering model with its restricted parameters as described above, we fitted 
the six available spectra. The results are presented in Fig.|9]in the form of unfolded spectra. It is 
worth warning that unfolded spectra should not be used to confirm spectral features in the data, but 
rather to present what a best- fit model spectrum would look like assuming it well represents the 
source spectrum. In each panel of Fig. |9l the diff"erent model components are plotted, namely the 
absorbed power law that turns over above 2 keV, the unabsorbed power law that turns over only 
below 1 keV due to Galactic absorption, and the reflection component. The relative contributions 
of the absorbed and unabsorbed components, i.e. the partial covering, varies dramatically between 
epochs. The absorbed component is very weak in the 2007 XMM-Newton observations (top panels 
of Fig. O, as the fit provides only upper limits to its presence and indicates very low covering 
fraction (< 20%), if any. Such an absorber would be consistent with the highly absorbed hard 
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X-ray component detected with Suzaku by iReeves et alj (|2009|) . but we can not conclusively 



confirm its presence in the 2007 XMM-Newton spectra. The reflection contribution to the total flux 
is generally small, but persistent in all epochs, and its emission below 1 keV is required to produce 
the soft excess observed in some of the spectra as well as the Fe-K line emission and the spectral 
drop off" above 5 keV. It can be seen from Fig. |9] that the models provide reasonable representations 
of the wide variety of spectra observed from PDS 456. Despite the several fixed parameters and 
most notably the power law slope, by varying the column densities and the fractional covering by 
the outflow, the quality of the fits is reasonable with;^f^ / d.o.f ~ 1.0 - 1.3, although not always 
statistically acceptable. We believe the imperfect fits are not necessarily due to the inadequacy 
of partial covering for describing PDS 456 at all epochs. More likely, it is the shortcomings of 
the models and their incomplete detail. In some cases, calibration issues can also afl^ect the fits. 
For example, ignoring the instrumental Si edge region of the Suzaku spectrum between 1.75 and 
1.9 keV dramatically improves / d.o.f from 1.30 to 1.20. 

The best- fit model parameters for the six spectra are summarized in Table [31 In light of 
the excessive local column density found in the 2007 XMM-Newton observations, the absorption 
in the model is allowed to exceed the nominal Galactic value of Nh = 2 x 10^^ cm"^. 
However, the slightly higher local column density detected in the 2007 XMM-Newton spectra 
(Nh = 2.7 X 10^' cm~^, see Section [3]) is not required in other epochs, except perhaps in the 
lower-quality ASCA spectrum where it is not as well constrained. As discussed in Section |3l this 
excess must be in our own galaxy (or at z = 0), but why it varies remains unclear. In the 2007 
XMM-Newton observations (top two panels in Fig.|9l), the unabsorbed component dominates (i.e., 
low covering fraction) and the outflow, if at all present, covers not more than 20% of the source; In 
the Suzaku, ASCA, and Chandra spectra, it is the absorbed continuum that dominates (i.e., almost 
full coverage), while the model for the 2001 XMM-Newton observation suggests an intermediate 
case in which about half of the source is covered by the absorber. Thus, the variation of the partial 
covering fraction can explain much of the spectral variability observed in PDS 456. This would 
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place the absorber at a distance from the source that is comparable to the source size. Variable 
partial covering, however, can also be due to many small sources, e.g., an assortment of magnetic 
flares above an accretion disk, some of which are absorbed and some of which are not. The size 
of the patchy absorber then would be of the accretion disk scale. 

Apparently, the column density in the absorption components also varies dramatically 
between epochs. A column of a few times 10^^ cm"^ or higher can account for the spectral 
curvature between 2-5 keV, which is observed in the 2001 XMM-Newton, the ASCA, and the 



Chandra spectra. The Suzaku spectrum requires only the tw o high-ionization ( 



og^ > 3.5) 



20091) and that cover 



absorbers that produce the blueshifted Fe-K absorption lines (|Reeves et al. 
most of the source. In the 2007 XMM-Newton spectra, the alleged absorber does not significantly 
affect the spectrum up to 12 keV. Indeed, the fits with absorption do not represent an improvement 
over those without absorption presented in Section [31 despite the high column densities (but low 
covering fraction). The absorption model in the 2001 XMM-Newton observation requires (at least) 
two components, and it needs to describe the absorption detected in the (simultaneous) RGS 
spectrum of the same observation. The most conspicuous absorption feature in that RGS spectrum 
occurs at approximately 0.9 keV and can be discerned with the CCDs. Such deep absorption lines 
require an absorber that essentially fully covers the source. On the other hand, partial covering is 
required to produce the observe d broad-band curvature of the spectrum. This result was already 



reported by 



Reeves et al 



(|2003|) . Hence, the model proposed here for the 2001 XMM-Newton 
observation, which is shown in Fig.|9l has a high-ionization absorber fully covering the source and 
a low ionization component only partially covering it (« 50%), both with columns of a few times 
10^^ cm~^. An additional component may be ne eded to fully accoun t for the deep and blueshifted 



Fe-K edge observed at 7 - 8 keV (see Fig. |9] and 



Reeves et al. 



20031) 



The reflection component is consistently present in all epochs. Its relatively moderate flux 
indicates the reflector only partially intercepts the X-rays from the primary power law source. In 
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fact, the models for all six epochs indicate it reflects about 20% of the primary source (see also 
Section [3]). This result is generally consistent with reflection off" a moderately coUimated outflow. 
However, if the outflow is patchy as suggested by the varying partial covering, the outflow 
segments can in fact be spread over a broader solid angle. In the 2007 XMM-Newton observations, 
the reflecting matter appears highly ionized (log ^ up to 4) and thus does not produce the soft X-ray 
emission that appears in all of the other epochs where it is approximately in the 2 < log^ < 3 
range. The reflection models tend to systematically prefer a super-solar Fe abundance, although 
the exact value is poorly constrained. We deem such a high Fe abundance to be unphysical, and 
ascribe it mostly to the shortcomings of the model, rather than to an actual abundance efi'ect. 
Since in any case, without multiple identifiable narrow features it is impossible to expect a reliable 
assessment of abundances, we fix the Fe abundance at three times solar (Section [3]). The outflow 
velocities of the absorber/reflector are mostly constrained by the soft X-ray emission lines in 
the reflection component, which yields effective redshifts in most cases of z = 0.10 - 0.14. If 
associated with an outflow in the reference frame of PDS 456 (z = 0.184), these redshifts translate 
into outflow velocities of Vom/c ~ -0.04 to -0.08. Two exceptions to this velocity range are noted. 
The first is the Suzaku spectrum in which t he narrow absorption lines require a much faster outflow 



of i^out/c = -0.26 and -0.31 as reported by 



Reeves et al. 



(120091) . The second exception is the 2001 



XMM-Newton observation that has a deep blueshifted edge-like feature at approxim ately 8 keV. 



When ascribed to the outflow, this edge implies Vout/c ~ -0.19 ([Reeves et al 



I 



2003 ). 



Although the reflection component of the model does produce the observed flux below 1 keV, 
it does not always match all of the spectral details very well. Consequently, at least formally, the 
fits are not always statistically acceptable. Indeed, invoking a partial covering absorber inserts 
additional uncertainty, as the fit is only partially sensitive to adding continuum and absorbing it. 
Nonetheless, we suspect the deficient fits are likely due to shortcomings of the detailed models we 
use, and do not necessarily invalidate the relevance of partial covering or reflection. For example, 
even moderate broadening of discrete reflection features, due to velocity or to scattering, would 
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make them less prominent. Also, incomplete atomic data, the constant-density approximation, 
and the simplistic geometry used in the reflection model certainly affect the goodness of the fit. 
Discrete features detected in the data, e.g. in the XMM-Newton 2001 and Suzaku 2007 spectra, raise 
the confidence one has in the model, even if they do not result in superior statistical fits. In fact, 
the 2001 XMM-Newton spectrum is perhaps the hardest to fit and we find the partial covering model 
essential to even obtain an adequate description of the continuum. In the following section we 
elaborate on the model for this observation that is derived from both the EPIC PN and RGS data. 



4.3. The 2001 RGS spectrum revisited 



The 2001 XMM-Newton observation is the only observation of PDS 456 to provide identifiable 
absorption lines in the gratings. The least ambiguous absorption features in the 2001 RGS 
spectrum are observed in the 1 1.5 - 14.0 A range, and in particular there appears to be a deep 
and broad absorption trough between 13.5 - 14.0 A. We re-analyze the RGS spectrum of PDS 456 
from the 2001 observation with the same absorptio n models as used ab ove for the CCD data. The 



Reeves et al. 



(12003b, who used XSTAR tab 



Reeves et al. 



(2003) 



es. 



method used here is not very different from that of |j 
We find that the ionization parameter of log^ = 2.6 (cgs units) given by 
for the RGS spectra produces on top of the 14 A feature other strong lines that are not seen in 
the data. In order to produce only a few features, we need to increase the ionization parameter 
to log^ Si 3.0, to ensure that most species are fully ionized and only a trace of Fe-L remains. 
Indeed, we are able to find a model with log^ = 3.1 that fits the RGS absorption features, and is 
marginally consistent with the fit to the EPIC PN data obtained in Section [42] (Table [3]) that is not 
as sensitive to narrow line features. The relevant part of the RGS spectrum along with the model is 
presented in Fig.[lOl Unfortunately, RGSl does not provide data between 1 1.5 - 13.75 A, due to a 
failed chip. In Fig. [lUl thus, only RGS2 is shown, although both RGSs are used in the fit and are 
consistent with each other. 
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According to the present model, several absorption features are shifted by z = 0.13 and 
can be tentatively identified in Fig. |9l We identify the wide trough between 13.5 and 14.0 A 
predominantly with Ne^^ Ly a (A^est = 12.134 A) that is broad (vtmb ~ 2500 km s"') and blended 
with the resonant (/ = 1.24) 2p - 3d line of Fe+^o (Arest = 12.284 A). The strength of the 
Ne^*^ Ly a line may indicate an overabundance of Ne, although no abundances could be reliably 
constrained. The Ne^^ Ly/3 line in this model falls at 11.6 A (A^c&t = 10.239 A), just next to 
a chip gap in RGS2 (and where RGSl is still inactive), so it is not obvious in the spectrum. 
Another relatively strong line is found at 12.01 A and identified as the Fe"^^^ 2s - 3p doublet 
(/Irest = 10.621, 10.664 A). Finally, the model plotted in Fig.[lO]also predicts absorption by Fe"^^^ 
at 12.4 A (/Irest = 10.982, 1 1.021 A), which is not clearly identified in the spectrum, but consistent 
with the data uncertainty. This blend could be overestimated in the model if the ionization 
parameter or the ionization balance are slightly ofl". These high charge states of Fe give an idea of 
the highly ionized state of the absorber. 

The outflow velocity resulting from the aforementioned tentative line identifications as due 
to a shift of z = 0.13 is Vn„t/c = 0.13-0.184, or Vout = -16, 000 + 1, 500 km s~\ slower than 



found by 



Reeves et al. 



(l2003h . but consistent with the velocities seen for HLya absorption in 



the UV. The connect ion to the UV absorb er is puzzling since no high-ionization absorption lines 



were detected there (lO'Brien et al 



20051). The determination of the outflow velocity obviously 



depends on the correct identification of the absorption lines. Except for the 14 A feature, other 
identifications are marginally significant, and thus the outflow velocity derived mostly from one 
feature will always remain ambiguous. The Fe-K edge observed with EPIC PN at 8 keV still 
requires a high velo city of ~ 50,000 kms This leaves the present velocity as well as that of 



Reeves et al. 



(|2003|) possible. In any event, the CCD spectra that are of higher statistical quality 
and that were discussed in the previous section can not distinguish between the two interpretations 
of the RGS data. Fig. [10] also highlights emission lines produced by the reflection model, which 



appear to be consistent with the RGS data. The brightest line is at 17 A, which is the 2p - 3d 
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resonance line of Fe"^ (redshifted with the outflow by z = 0.13). Two additional narrow features 
that are not reproduced by the model appear in the RGS2 spectrum of Fig.[10l namely at 15.18 A 
and at 17.90 A. The 15.18 A line appears in both RGSs and corresponds to 13.43 A in the outflow 
frame, which could be the Is - 2p resonant line of Ne^^. The 17.9 A line is totally absent from the 
equally good spectrum of RGS 1 (not shown) and can be discarded. 

5. Conclusions 

This paper presents the XMM-Newton observation of PDS 456, which took place in September 
2007, and attempts to generalize its results to the broader context of four other CCD-quality 
observations of the quasar acquired over the past decade, each one of them revealing a spectrum 
quite remarkably different from the others. In a sense, the 2007 XMM-Newton spectrum is the most 
straightforward to interpret, as no unambiguous absorption is present. On the other hand, owing 
to the lack of ionized absorption, a significant reflection component can be discerned. Global fits 
to the spectrum suggest it may be outflowing with respect to the quasar's rest frame, although this 
result is less robust that the mere existence of a reflecting medium. Since an absorbing outflow has 
been identified in PDS 456, a natural conjecture is to associate the reflector with outflow segments 
outside the line of sight. The small solid angle coverage of the source of ~ 20% found for the 
reflecting medium is consistent with this notion, although the assumption of a single scatter, 
perfect reflector are clearly an over simplified picture of a more realistic reflecting outflow. The 
velocities required to model the reflecting material are in the range of v/c = -0.04 to -0.08, which 
is roughly comparable to the blueshifted absorption measured in the UV and to the interpretation 
of the 2001 RGS spectrum. A separate component is needed to explain the much faster Fe-K 
absorption, observed e.g. with Suzaku. 

The historical spectra of PDS 456 point to a complex nature of the source and complex 
absorption. We tested whether all spectra can be explained by a partial covering outflow model 
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that both absorbs and reflects the central source and obtained mixed results. This type of model 
is sufficiently flexible to produce a wide variety of spectra, but is only partially successful in 
reproducing those of PDS 456. While the diff'erent over all spectral curvatures of PDS 456 are 
produced by varying the fractional coverage of the source essentially from zero to unity, the 
spectral fits are not always statistically acceptable, and some narrow features predicted by the 
models are discrepant with the observed spectra. These features vary in the data on time scales 
from days to years. The reflection component can explain some of the emission that is excessive 
of the power law continuum. In the 2007 epochs, it provides a good fit to the Fe-K emission, 
but in other epochs, its details of the soft emission below 2 keV do not always match the data 
very well. Future X-ray instruments that would provide high S/N spectra on short time scales are 
imperative to further advance the study of such outflows. 

EB acknowledges funding from a NASA XMM-Newton Guest Observer grant and from NASA 
grant 08-ADP08-0076. SK is supported at the Technion by the Kitzman Fellowship and by a grant 
from the Israel-Niedersachsen collaboration program. 
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Fig. 1. — EPIC PN light curve of PDS 456 from two consecutive XMM-Newton orbits during 2007, 
Sep. 12-15. Data are binned to 500 s and include the 0.2-15 keV photon energy range. 
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Fig. 2. — Typical ofF-source EPIC PN background spectrum from the XMM-Newton observation of 
PDS 456 during orbit 1421. 
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Fig. 3. — Ratio of PDS 456 spectra as observed by the EPIC FN camera in Sep 2007 during XMM- 
Newton orbits 1421 (black) and 1422 (red) to a Galactically absorbed {Nh = 2.0 x 10^' cm"^) 
power law (top), with an increased Galactic column density of Nh = 2.7 x 10^' cm~^ (middle), 
and after adding a reflection component (bottom). All spectra are binned to 40cr significance, but 
not more than 20 bins, and are plotted on the same scale to enable easy comparison. The resulting 
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Fig. 4. — {top) EPIC PN spectra and best-fit absorbed power law models (c.f., Fig. [3] middle panel). 
{bottom) Contributions to from XMM-Newton orbits 1421 and 1422 around the Fe-K region. 
Spectra are binned to AOcr significance, but not more than 20 bins. Note the overall similar excess 
emission from ~ 5 to 7 keV and the difference in some spectral details between orbits. 
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Fig. 5. — Two-parameter statistical confidence contours of the ionization parameter (^) and redshift 
(z) of the reflection component fitted to the 2007 EPIC PN spectra of PDS 456. The three contours 
represent b^x^ of 2.3, 4.61 and 9.21 with respect to the best values of ^ = 10000 and z = 0.122. 
Values of ^ > 10000 are not available with the current reflection model, but the high confidence in 
the result of z < 0.184 (i.e. outflowing reflector) appears to be robust. 



-31 - 



<M - 



CO 
DC 



CD 
T3 
O 



CO 

-I— » 

CO 
Q 



<M - 



<M - 



C\J - 




iiikijliLAiiij.j4i..iiiliiliLJiilrtiiL. 



10 



15 20 
Observed Wavelength (A) 



25 



Fig. 6. — Ratio of PDS 456 RGS spectra to the best (continuum) model obtained from fitting the 
EPIC PN spectra. Data are shown separately from top to bottom for XMM-Newton orbit 142 1 : RGS 1 , 
RGS2, and for orbit 1422: RGS 1 , and RGS2, and are rebinned uniformly by a factor of two. Straight 
line segments with no errors are due to gaps in the RGS detectors. No absorption (i.e., outflow) 
features recur consistently in several data sets. 



-32- 




Observed Wavelength (A) 

Fig. 7. — Zoom into ratio of PDS 456 RGS spectra to the best (continuum) model obtained from 
fitting the EPIC FN spectra. Data are shown separately from top to bottom for XMM-Newton orbit 
1421: RGSl, RGS2, and for orbit 1422: RGSl, and RGS2, and are not rebinned beyond the original 
reduction (Section [21). Straight line segments with no errors are due to gaps in the RGS detectors. 
No compelling evidence for recurring absorption (i.e., outflow) is seen. 
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Fig. 8.— PDS 456 as observed by EPIC PN in Sep. 2007 (black), Suzaku in 2007 (red), EPIC PN 
in 2001 (green), Chandra/MEG in 2003 (blue), ASCA/SIS in 1998 (magenta), and RXTE/PCA in 
2001 (cyan). The ratio of data to a Galactically absorbed power-law spectrum with photon index 
r = 2 is plotted for all observations. 
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Fig. 9. — Unfolded spectra of PDS 456 sampling its spectral variability over the past decade (in 
arbitrary order). Spectra are produced using the partial covering model discussed in the text. The 
three dotted lines in each spectrum delineate the absorbed, and unabsorbed power law components, 
and the reflection component. The resulting parameters are summarized in Table [3l 
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Fig. 10. — Extract forni RGS2 spectrum of PDS456 obtained in 2001, and the proposed XSTAR 
absorption model (solid line). This model corresponds to z = 0.12 (an outflow velocity of 16,000 
km s~^ in the quasar frame) that produces the most significant (redshifted) absorption feature at 
13.7 A as a blend of Ne""^ Ly a and Fe^^" 2p - 3d, the line at 12 A as the Fe"'^^ 2s - 3p doublet, and 
the 2p - 3d emission line of Fe"^^^ at 17 A from reflection. 
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Table 1. Observation Log (in arbitrary order) 



Telescope Instrument Obs ID. Start Date Exposure (ks) 



XMM-Newton 


EPIC PN 


0501580101 (orbit 1421) 


2007 Sep. 12 


91 


XMM-Newton 


EPIC PN 


0501580201 (orbit 1422) 


2007 Sep. 14 


88 


XMM-Newton 


EPIC PN 


0041160101 (orbit 223) 


2001 Feb. 26 


44 


Suzaku 


XIS 0, 3 


701056010 


2007 Feb. 24 


180 


Chandra 


HETGS 


4063 


2003 May 7 


143 


ASCA 


SIS 


76078000 


1998 Mar. 7 


43 
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Table 2. Model parameters obtained from fitting 2007 EPIC PN spectra 

Parameter Orbit 1421 Orbit 1422 

local (10^2 cni-2) 
POWER LAW 
photon slope 

flux density at 1 keV (10"^ s'^cm-^keV"!) 
IONIZED REFLECTOR 
Fe abundance (solar units) 
slope 
redshift 

ionization parameter ^ (erg s"^cm) 
angular coverage 
d.o.f. 

Ar'/d.o.f. 

^Only one value is quoted for parameters that were tied for both data sets. 
Errors are 90% confidence limits. 

''Redshifted uniformly to the PDS 456 rest frame of z = 0. 1 84. 

"^Values greater than 10000 erg s"^cm are not available in the model. 

'^Assuming no losses on the reflector. 



261^^-^^ 

"•^"'-0.002 



2.20 ± 0.01 2.25 ± 0.01 

4 0+0.2 2 4+0.2 



3 3 "^'-^ 

-0.6 

tied to power-law slope 
12 

"•^^ -0.01 

10000 

0.19±0.07 0.27±0.10 
2356 
1.05 
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Table 3. Best-fit parameters for partially covered absorption models with reflection 



Parameter XMM-Newton XMM-Newton XMM-Newton Suzaku ASCA Chandra 

EPICPN EPICPN EPICPN XIS0&3 SIS HETGS 

2007a 2007b 2001 2007 1998 2003 



local (10^2 cm-2) 




26"'°°2 


21+0"2 


21+''"2 


3+°-i 


20'^°-i 


POWER LAW ^ 














photon slope 






fixed at 2.25 






1 keV flux density 


5.1 ±0.3 


2.5 ±0.1 


6.2 ± 0.2 


3.1 ±0.3 


4.6 ± 0.9 


3.2 ± 0.2 


(10-3 s-icm-2keV-i) 














IONIZED REFLECTOR 














photon slope 






fixed to the primary power law 






^ (erg s-'cm) 


10000 

-2000 


10000 +" 

-4000 


1250 

i^J^ -100 


1200 


850 +300 

-150 


400 +200 

-100 


Ape (solar units) 






fixed at 3.0 








ang. coverage Q,-efl/27r 


0.12 + 0.07 


0.26 + 0.10 


0.14 + 0.05 


0.20 + 0.07 


0.08 + 0.08 


0.12 + 0.12 


IONIZED ABSORBER 














log^high 


3.6 




3.0 + 0.1 


> 3.5 


2.7 


1.5 


<'^\l022cm-2) 


100 




2 ±0.5 


100 + 30 


100 ± 30 


0.2 ± 0.2 


log^i"™ 


-2.7 


-1.7 


-1.5 


> 3.5 


-3.2 


-1.8 




100 


100 


6±2 


100 ± 30 


3 ±3 


8 ±4 


covering fraction 


<0.2 


< 0.1 


0.5 ± 0.05 


0.7 ± 0.3 


0.9 ±0.1 


0.75 ± 0.05 


z (refl/abs) 


0.10 ±0.02 


0.12 ±0.02 


0.13 ± 0.02 " - 


-0.08 ± 0.02 f 


0.11 ±0.02 


0.14 ±0.02 


Vout/c = 2-0.184" 


-0.08 ± 0.02 


-0.06 ± 0.02 


-0.05 ± 0.02 ^ - 


-0.30 ± 0.03 f 


-0.07 ± 0.02 


-0.04 ± 0.02 


d.o.f 


1281 


1066 


1053 


214 


105 


275 




1.10 


0.99 


1.18 


1.30 


1.00 


1.29 


0.4-12 keV model flux 


10.9 


6.7 


11.7 


7.8 


4.7 


4.8 



(10-12 erg s-icm-2) 
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^Shifted to quasar rest frame at z = 0.184 
''Assuming no losses on the reflector. 

"^Fixed absorber properties: Vturb = 2, 500 km s"^ and solar abundances. Values without errors are not constrained. 
''Conamon to reflector and absorber, except for Suzaku where absorber velocity is well constrained from Unes. 
^For reflector and fuUy covering absorber 
^Independent redshift for the reflector of 0.06, or Vout/c = -0.12. 



